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q\ ■ Abstract 

We show that the lowest-order QCD calculation in a simple model of 
o ■ elastic vector-meson production does reproduce correctly the ratios 

of cross sections for p, <fi and J/ip, both in photoproduction and in 
high-Q 2 quasi-elastic scattering. The dependence of the slopes on 
the mass of the vector meson is reproduced as well. We examine the 
lower-energy data, and find that the energy dependence of the cross 
section does not depend on Q 2 , but may depend on my. 
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Elastic vector-meson production opens a precious window on the interface 
between perturbative QCD (pQCD) and non-perturbative hadronic physics, and 
is complementary to deep-inelastic scattering (DIS). Indeed, DIS has been until 
recently the triumph of perturbative ideas, whereas elastic processes were mostly 
treated through non-perturbative methods. These two processes now meet at 
HERA, where pQCD and Regge models have to be merged to obtain a full un- 
derstanding of the data. Elastic vector-meson production has the extra advantage 
of containing by definition two scales, the mass of the produced vector meson, 
My, and the off-shellness of the photon, Q 2 = —q 2 . We shall see that both 
dependences can be understood through a lowest-order calculation. 

The data for j*P ~* V"Pi which both HI [I], || and ZEUS have obtained, for 
V = p, (ft and J /if), exhibit the following main features: (i) the Q 2 -distribution of 
the cross section is shallower for J / if) than for p, the two cross sections becoming 
comparable around Q 2 = 10 GeV 2 and (ii) the t-slopes of the differential elastic 
cross sections depend on the mass of the produced meson, and become shallower 
as the mass increases. 

Several models have been proposed to describe this process. Originally, 
Donnachie and Landshoff (DL) Q extended their soft-pomeron model to predict 
p production at EMC ||. This model works well there, but it has not been 
applied to predict the ratios of produced vector mesons. In photoproduction, 
DL preferred to resort to the quark counting rule to predict the ratios of cross 
sections @, and reached the conclusion that it does not work perfectly, but argued 
that the violations were reasonable. The transition to perturbative QCD was first 
introduced by DL 0], who noticed the analogy between the pomeron expressions 
and two-gluon exchange, in the transverse case. They used a "constituent gluon" 
propagator and two-gluon exchange to model the pomeron. This analogy was 
pursued by one of us f|, who showed that such a model can give reasonable 
agreement with EMC data. This was later confirmed by NMC |J. The final step 
to pQCD was performed by Ryskin |J who observed that at high-Q 2 and high 
My, the effective intercept should be analogous to that found in xg(x). This was 
later confirmed by Brodsky et al. fli~0" |. So far, there has not been a model applied 
to the full range of masses and Q 2 : this is the object of this letter. 

A model for exclusive vector-meson production must include three sub-models, 
for which we shall adopt the simplest ones: the transition 7* — > V is described 
by a zero Fermi momentum wavefunction, the colour-singlet exchange is modeled 



a la Low-Nussinov (TT|, and we shall only consider the constituent quarks of 
the proton. There are in principle 72 diagrams contributing to the amplitude: 
the gluons can be hooked 4 different ways to the quarks of the vector meson, 9 
different ways to the quarks of the proton, and both the direct and the s <-» u 
channels contribute to the amplitude in the high-s limit. As we shall explain 
however, the calculation of each part of the amplitude can be greatly simplified, 
so that one needs to calculate only the two diagrams shown in Fig. 1. 

Let us first consider the kinematics of the process, which is spelled out in 
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Figure 1: The two diagrams accounting for the transition 7* — > V. The 
dashed line represents the cut which puts the intermediate state on-shell. 



Fig. 1. Let P and q be respectively the 4-momenta of the proton and the photon, 
and P — A and q + A the momenta of the final-state proton and vector meson. 
We work in the high-w 2 limit, with w 2 = (P + q) 2 . The on-shell condition for the 
proton and the vector meson imply that A is transverse to order 1/w 2 : A Pt Ay 
with A r • P = A T ■ q = 0, and |A|| « |A 2 | = -i. 

As it will turn out, the imaginary part of the amplitude is proportional to 
w 2 . Crossing symmetry and analyticity then imply that the amplitude is purely 
imaginary, up to terms of order 1/w 2 , because the exchange is C = +1. Hence 
in the following, we shall calculate only the imaginary part, using Cutkovsky's 
rules, and putting intermediate states on-shell, as shown in Fig. 1. The quarks 
that make the vector meson are in the direction of q, whereas those that make the 
proton are parallel to P. The intermediate states come from the absorption of 
one gluon. In order for these states to be on-shell, the gluon must have vanishing 
components both in the direction of P and q, hence the gluon momenta are 
transverse to order 1/w 2 : k pa kr, with kx-P = kx-q = 0. 

Hence the gluon momentum is essentially transverse, and furthermore we do 
not need to worry about the (purely real) crossed diagrams. The transversality 
of the gluons implies that the embedding of a process at the quark level into a 



proton is particularly simple. We follow here Gunion and Soper |12| , and choose 
to represent the proton by a constituent model. This naturally leads to the quark 
counting rule, and has the extra advantage that the form factor of the proton is to 
a large extent measured. Indeed, when both gluons hit the same quark line, the 
form factor is given by the Dirac elastic form factor, measured in ep scattering: 

5 (f _ a2\ ^ (3-53 - 2.79t) 

If the gluons hit different quark lines, then the form factor is not known, but we 
know that infinite wavelength gluons must average out the colour of the proton, 
hence we know that this form factor has to cancel the infrared singularity that 
would result from the pole in the gluon propagator, and therefore it must reduce 
to £\{t) when either gluon goes on-shell. We choose to parametrise this form 



factor as: 

£ 2 (k, k — A) = E x {k 2 + (k - A) 2 + c k ■ (k - A)) (2) 



There are theoretical arguments |13| as well as phenomenological ones [H which 
lead to the conclusion that c ~ — 1. 

The rule is then to calculate the process at the quark level. This leads to 
an integral over the transverse momenta of the gluons. One then introduces the 
difference of form factors 3(£i — £2) into the integral to get the same process at 
the proton level, thereby reducing the number of diagrams by a factor 9. For 
the lower trace, and in the high-?/; 2 limit, we need to keep only the leading terms 
in p, the momentum of the quark, hence the trace along the quark line in the 
square of the amplitude, is given by (4p a p f3 ) x (4p a p 13 ) (including a factor of 
1/2 for spin averaging), with p the momentum of the quark inside the proton, 
and and f3^'' the indices of the gluons. Thus we can treat the process at the 
level of the amplitude, without the need to square it, provided that we write the 
contribution of the lower quark line as 4p a p^. 

For the vector meson, we use a different model than for the proton, as we want 



to take into account the mass of the meson. The vector meson is modeled |15 
by its lowest Fock state, with no Fermi momentum, which implies m q = my/2, 
and the Vqq vertex, including the two quark propagators, is given by: 

$y = C (m q - 7-e (m q + 7.1; ) = C 7-e (m v + 7.V) (3) 



where v=^ is the quark momentum within the vector meson and C = J /y/24 is 
the normalisation that reproduces the vector meson decay rate, with 
fv ~ 0.025(GeV) m v the vector-meson decay constant squared. This effective 
vertex includes the propagators of the quark lines flowing into it. Reversing the 
direction of the quark current gives the same contribution, therefore we end up 
with only 2 diagrams to calculate - those shown in Fig. 1. 

The traces corresponding to the upper bubble, dotted into p a P/3, are: 

Ti = Tr[<£ (7.1/ + m) j> (7.(17 - v + k) + m q ) p 1 (7.(5 - v) + m q ) f\ (4) 
T 2 = Tr[^ (7.V + m) j> (7.(5 — v — k + A) + m q ) ( 

x(7.(-v - k + A) +m q ) f[ (5) 

with e and e the polarisations of the vector meson and of the photon. The 
propagators of the off-shell quarks are: p\ = —{Q 2 + my — t)/2 and p 2 = (rriy + 
Q 2 — t — Ak.q)/2. The answer is then proportional to: T = Ty/pi + Tijvi- This 
answer is explicitly gauge invariant: substituting e = q in T gives 0[[ 

In the transverse case, the leading dot product is simply w 2 = (p + q) 2 ~ 2p.q. 
For longitudinal polarisation, further contributions appear: e L .p w — tD 2 /(2v / Q 2 ") 

^Note that this was not the case for the expressions previously published in ||. The difference 
at high-Q amounts to a factor of 2 in the amplitude, presumably the discrepancy pointed out 
in M. 
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and cl-P ~ — w 2 /(2mv)- Keeping track of these, the leading term is proportional 
to: 

_ 2m v k ■ (k - A) 



t - ml - Q 2 + 4 k ■ (k - A) 

/ 2A.e e.p + 2g.e e.p — e.e w) 2 2 \ 

x ( (■„;.. + # e - pe - p J (6) 

Putting everything together, we obtain the following expression for the am- 
plitude: 



A =- ia stf m l\j^J ¥ 0^- 2 Sm)-S 2 (k,k-A)]x32T (7) 

where gy m = £,y/4ira e i m is the electromagnetic coupling of the different vector 
mesons: £ = 4g for the p, —1/3 for the and 2/3 for the J/^- For the various 
possible helicities, Eq. (^) gives: 

*A(T — >■ T) = ^2644#mvvW5 

/• d 2 fc r [g 1 (t)-g 2 (fc,fc-A)] [fc-(fc-A)] 

7 P(A; - A) 2 (t - m 2 v - Q 2 + 4 k ■ (k - A)) (m 2 v + Q 2 -t) 1 J 

A(L — >■ L) = ^x^(T^T) (9) 
my 

The helicity violating amplitude ^4(L — > T) is suppressed by 1/w 2 . As previously 
advertised, this answer is proportional to w 2 , therefore to w 2 , and is thus purely 
imaginary. The resulting cross section is independent of w 2 . Clearly this model 
cannot say anything about the energy dependence of the cross section. We shall 
assume that it comes in as a factor, R, and check whether the latter is mass- or 
Q 2 -dependent. Note that we can only determine the value of that factor times 
a 2 s . In the following, we shall let as run with the off-shellness of the gluons, 



and freeze it at some value a s . However, as we shall see later, the dominant 
contribution comes from gluons of small off-shellness, and the results we obtain 
are identical to fixed-coupling results for as = a^- 
The differential cross section is given by: 



da (Lot d<JL R 



+ e 



dt dt dt 16tt(w 2 ) 2 



\A(T — T)\ 2 + e\A{L -> L)\ 2 (10) 



with e the polarisation of the photon beam: e ~ 1 at HERA and e ~ 0.75 at 
NMC. 

We first give the results that we obtain for the various cross sections measured 
by ZEUS and HI. First of all, we show in Fig. 2 the dependence on Q 2 and my of 
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Figure 2: (a) Cross sections as functions of Q 2 , compared with data from 
HI 0, |] and ZEUS g at < w >«100 GeV, (b) Ratio of cross sections as 
functions of Q 2 at < w >~100 GeV, compared with data from HI ^] 
and Zeus HI. 



the integrated elastic cross section a(Q 2 ). We see that a common (Regge) factor 
is consistent with the data taken at HERA, as shown in Fig. 2(a). We insist 
on the fact that this factor is independent both of Q 2 and of my, as one would 
expect within Regge theory. Selecting only high-Q 2 data leads to a best value 
Ra 2 s = 0.6, with a x 2 /d.o.f. = 0.39. We do not find that Ra 2 s varies significantly 
within the energy range of HERA. Although we see no reason why our model 
should work in photoproduction, it turns out that our curves do go through the 
photoproduction points. Including these in the fit brings Ra 2 s to 0.55. 



meson 


A (nb) 


n 


P 


5534 


2.5 


<P 


1035.9 


2.4 




269.53 


1.5 



Table 1: Result of a fit of a(Q 2 ) to A{Q 2 )~ n , for 5 GeV 2 < Q 2 < 25 GeV 2 . 

We give in Table 1 the result of a fit to a power of Q 2 at large Q 2 . Although 
asymptotically all cross sections behave like {Q 2 )~ 3 , in agreement with 



TU| , we 

see that the data collected at HERA are not yet in that asymptotic regime. Note 
that our calculation holds only for Q 2 << w 2 , and hence it is not clear whether 
the asymptotic regime will be reachable at HERA. Our model is equivalent to 
those of refs || with a gluon distribution corresponding to 2-gluon exchange 
g(x) ~ 1/x. We see that here the Q 2 dependence is already reproduced, hence 
there is no room for an extra dependence coming from the gluon distribution, 
and the use of the asymptotic formula [I7| is misleading. 
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In Fig. 2(a), we show both the systematic and the statistical errors. In the 
ratio of cross sections, some of the systematic uncertainties cancel, and the re- 
production of that ratio is a more stringent test of our model, especially as the 
normalisation then drops out of our prediction. We show in Fig. 2(b) the re- 
sult of such a comparison. Again, we see that our model fares well, even in 
photoproduction. 

Hence we can understand both the my- and the Q 2 -dependence of the cross 
sections. One might object that this is because these are concentrated at low t, 
and argue that the t-dependence has to be wrong, as this is one of the well-known 
problems of perturbative calculations applied to diffractive scattering. 

The behaviour of the slopes as a function of my and Q 2 can be understood as 
follows. We can approximate the proton form factor £\{t) — £2(^5 k — A) as being 
proportional to £\{t) [k • (k — A)], using the fact that F\ is close to an exponential, 
and expanding for small k-(k-A). The amplitude of Eq. (]?]) can then be written 

as C(my,Q 2 )J R(t)£i(t)J- ( 2 * z J, with C a constant with respect to t, and 



T a calculable function. This means that the logarithmic derivative of da/dt 
becomes: 

, , dSi dR 1 dJ~ , , 

W = Sjdt + Rdt + [m\ + Q 2 ) T 2 dx *= ; * l ! 



v 



Thus the slope is approximatively made of three terms: one corresponding to the 
proton response, one to the pomeron response, and one to the response of the 
loop which converts the photon into a vector meson. Only the latter depends on 
my and Q 2 , and we see that it decreases rather fast with both of these factors. 
At large Q 2 + m V) it becomes negligible, and only the first two responses matter. 
This means that this kind of model predicts that all the slopes have to reach the 
same asymptotic value. This value is about 4 GeV~ 2 . 

This variation of the slopes with Q 2 + my enables us to reproduce the mea- 
surements of HERA: in photoproduction, we obtain 1/ < t p >= 8.75 GeV -2 , 
1/ < t<p >= 7.35 GeV -2 and 1/ < tj/^ >= 4.55 GeV -2 , whereas the p slope is al- 
ready down to 5 GeV -2 for Q 2 = 5 GeV 2 . This is inconsistent with the data from 
HI, but agrees with those from Zeus. We want to point out that the experimental 
evaluation of the slope demands that the cross sections be exponential in t. We 
have quoted our results for 1/ <t>, which in the case of an exponential fall-off 
e bt is equal to the logarithmic slope b. The fact that the differential cross section 
is not an exponential makes the comparison with data difficult. To illustrate 
the effect of the curvature, we compare in Fig. 3 our results with the data for p 
photoproduction in HI. We see that although the slope is supposed to be 10.9, 
our curve reproduces the data fairly well, with a 1/ < t > of 8.75 GeV~ 2 . Hence 
the different definition of the slope leads to a possible systematic correction of 
about 2 GeV~ 2 . No matter which model is used, da/dt is not an exponential, 
and we urge the experimentalists to quote a < t > instead of a logarithmic slope. 
In fact, our model works although it assumes no shrinkage. There is room in the 
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p case to add a slope coming from the Regge factor R, which would contribute 
about 4.6 GeV -2 to b for a pomeron slope a' = 0.25 GeV -2 . One would then 
expect a similar contribution in the J/ip case. Hence if the HI measurements of p 
slopes are correct, it is likely that the J /if) numbers have a background problem: 
inelastic contributions would significantly reduce the slope. 

The only problem at HERA seems to be the helicity structure of the cross 
section. The data support the prediction that helicity is conserved, but the ratio 
ul/^t does not follow the results of our model. This ratio has to behave as Q 2 
for near-shell photons, as a consequence of gauge invariance. Our model fulfills 
this requirement, but it also predicts that this linear behaviour continues for 
all Q 2 . The data, on the other hand, seem to indicate that the ratio reaches a 
plateau around 2 at high Q 2 . This would indicate that our high-Q 2 transverse 
cross section is wrong. This is indeed possible: we have assumed that the meson 
wave function is dominated at high Q 2 by configurations in which both quarks 
have equal momenta. In fact, it is likely that further configurations exist ||18| 



which would give additional contributions to the transverse cross section. This 
may account for the fact that this model does not reproduce the ration ul/ut 
measured at HERA. We plan to examine the role of Fermi momentum in a later 



paper. Another interesting possibility has recently been pointed out [17|: the 



form factor describing the recombination of quarks into a p meson is taken to be 
1 once the quarks are restricted to the mass region of the p, and it is argued that 
the process will remain largely elastic. The ratio <Jt/&l is then shown, in the 
leading-log approximation, to be distorted, at t — 0, by the x and Q 2 dependence 
of the gluon distribution, in such a manner that it reproduces the data. This 
is the only model so far which manages to reproduce the observed ratio, but it 
can be used only for p mesons at t = 0, and heavily relies on the perturbative 
evolution of the gluon distribution. 

Finally, we can now examine the w 2 -dependence of the cross sections. We 
have seen that at HERA, the Regge factor does not seem to depend either on the 
meson mass or on Q 2 . This is clearly reminiscent of the behaviour expected from 
a simple pole. We adopt the same philosophy when fitting to lower-energy cross 
sections from EMC and NMC [H], and write R = s 2a °~ 2 , «o being the intercept of 
the pomeron. We show in Fig. 4 the curves which correspond to «o = 1.16. This 
value of the soft intercept is rather high, but may not be entirely excluded by 
fits to total cross sections, especially once the effect of unitarisation is taken into 
account |T9| , [20|| . We see that the fit is reasonable for the p and the 0, but fails 
for the J/ip. The fact that the p data is somewhat high can be understood from 
the fact that there are contributions from lower trajectories to the p production 
cross section. The interference between a/f exchange and pomeron exchange 
could contribute as much as 20%, hence a high p cross section. In the J /if) 
case, one has to realise that the calculation we have presented here is valid if 
Q 2 + rriy << w 2 . This is not the case for the J/ip at EMC. Besides, the data for 
J / if) production from EMC is not subtracted for any diffract ive background, and 
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Figure 3: Photoproduction differential cross section da/dt, compared with 
HI data §]. 



can only constitute an upper bound on the elastic cross section. It is likely that 
there is a contamination of the order of a factor 2 from inelastic contributions, 
as the inelastic background is larger than in the p case, for which the EMC data 
were severely contaminated ||. 

To check whether the behaviour is indeed that of a simple pole, we can fit the 
factor R separately for each meson. We give in Table 2 the best factors, and the 
corresponding intercepts, given the value of Ra 2 s at HERA. We do not give the 
results for the J/ip as it is impossible to obtain a good fit to the data. One sees 
that the <fi and p data require somewhat different factors, This can be seen more 
clearly when one notices that the ratio r of cross sections has changed when going 
from HERA (r = 0.18) to NMC (r = 0.12). Hence there may be a problem with 
the smallness of the (f> cross section at NMC, which seems to imply an intercept 
too large to be compatible with that of a soft pomeron. On the other hand, there 
is no sign in either p or data of a Q 2 dependence of the intercept. 



meson 


Ra 2 s 


X 2 /d.o.f. 


energy (GeV) 


intercept 


P 


0.4 


0.6 


40 < v < 180 


1.10 





0.268 


0.138 


40 < v < 180 


1.206 


combined 


0.32 


1.19 


40 < v < 180 


1.161 



Table 2: best values of Ra 2 s at NMC. 



Before concluding, we must mention that although the above looks like a 
successful perturbative calculation, most of the contribution to the total cross 
section comes from the infrared region. Imposing a cut-off \k 2 \, \ (k — A) 2 | > 
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Figure 4: Cross sections compared to lower-energy EMC and NMC data |], 
at < w 2 >«200 GeV 2 . 



1 GeV 2 reduces the cross section by a factor 10. This dominance of the infrared 
region justifies a posteriori our choice of scale in as' we see no theoretical reason 
to make it run with either Q 2 or rriy, as these scales are unrelated to the off- 
shellnesses of the gluons entering the vertices. The dominance of the infrared 
region confirms the results of []2T] , where it is argued that the BFKL perturbative 
resummation is dominated by the non-perturbative region. 

The simplest modification to the infrared region follows the ideas of 
Landshoff and Nachtmann p2[ , which have recently been further motivated by 
lattice studies that the gluon propagator needs to be modified at low k 2 , 
taming its behaviour to something softer than a pole. One of the main effects of 
that modification is that the differential cross section becomes much more linear, 
and that as a consequence < t > becomes bigger. We have checked that such 
a model has all the other features of the one detailed above. The only major 
difference is that the slopes reach their asymptotic values much sooner, but it 
remain impossible to accommodate both the p shrinkage and the J/ip data. 

To sum up, we have shown that many features of elastic vector-meson pro- 
duction at HERA can be understood in a simple QCD model. We have also 
seen that comparison with lower-energy data proves somewhat problematic. The 
HERA data seem to indicate that the Regge factor does not depend either on Q 2 
or my. The comparison with NMC p data points to a soft pomeron intercept, 
whereas the <ft data would favor a larger one. We do not believe that much can 
be concluded from the J ftp data, which seem to have an inelastic background. 
Following [14] we plan to model this background in a further publication. 
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